THE POLYA-TCHEBOTAROV PROBLEM 
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Abstract. We describe the solutions to the problem of identifying the con- 
tinuum in the complex plane that minimizes the logarithmic capacity among 
all the continuum that contain a prefixed finite set of points. This description 
can be implemented numerically and this can be used to improve the estimates 
on the Bloch-Landau constant and other related problems as the maximal ex- 
pected lifetime of the Brownian motion on domains of inner radius one or the 
principal eigenvalue for the Laplace operator on such domains. 



1. Introduction and history of the problem 

Polya in [P6129j discussed the following problem which was suggested to him by 
Tchebotarov: 

Problem 1. Given a finite number of points E := {ai,...,a„} C C, find the 
continuum K with minimal logarithmic capacity such that E d K . 

For any continuum A', its complement in the Riemann Sphere, = Coo is simply 
connected, therefore there exists a unique conformal map f : fl ^ Coo \ D{0, R) 
such that f{z) = z + co + c_i/z + • • • around oo. Here R = R{n,oo) is called 
the conformal radius of O with respect to oo. In fact, cap(if) = R{n,oo). This 
provides an equivalent formulation of Problem [H which is usually called the outer 
reformulation of the Polya-Tchebotarov problem: 

Problem 2. Given a finite number of points E := {ai, . . . , a„} C C \ {0} find a 
conformal map / : D ^ C \ i? such that /(O) = and |/'(0)| is maximal. 

So we are looking for a simply connected domain Q that contains the origin, it is 
contained in C \ E' and such that the density of the hyperbolic metric at the origin 
is minimal. Such domain will be called an extremal domain and the corresponding 
conformal map, an extremal map. 

The existence of the solution is obvious by a normal family argument. This 
problem was studied in detail by Laurentiev in |Lau30| . He proved the uniqueness 
and the basic structure of the solution by the method of variations of the boundary. 
The structure of the extremal domain is characterized by the following theorem, 
see |Lau34| . 

Theorem 1 (Laurentiev). Given a finite number of points E = {ai,...,a„} C 

C, there exists a unique extremal domain VL = /(B) for the problem and it is 
characterized by the following properties: 

(1) Each point of the plane belongs to either Q, orV:— dVL. 
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(2) The boundary T consists of finitely many simple arcs of analytic curves. 
The points Oi and {00} are endpoints of n + 1 distinct arcs. Every point 
ofT different from the Oi or {00} either belongs to a unique arc and it is a 
regular point ofT, or it is the common end of at least three arcs. 

(3) To any arc afi consisting of regular points of T there correspond under the 
conformal mapping f~^ two arcs of the same length on the unit circle. 

When the last property [3] holds wc say that the arcs are harmonically symmetric 
with respect to the origin and it will be the key property to find a numerical 
algorithm to determine the solution to the problem mentioned above. 

In the proof of this last theorem, Laurentiev assumed that the desired domain 
is bounded by finitely many simple Jordan arcs. This assumption was removed by 
Goluzin who used the method of inner variations to prove the following: 

Theorem 2 (Goluzin, |Gol46| ). Let ai, . . . , a„ be arbitrary given points in C. Let 
K be the extremal continuum for Problem [7J Then K is the union of the closures 
of all critical trajectories of the quadratic differential 

Q{z)dz' = -flpii^dz^ 
U.k=iiz-ak) 

where bi are some unknown parameters. The extremal univalent function g : Coo \ 
D Coo \ {o-i, . ■ . ,o,n}, with 5(00) = 00 that maximizes g'{oo) must satisfy the 
following differential equation 



Remark 1. The points bi correspond to common end points of several arcs. If some 
point hi is a common end of m arcs, then the term g{z) — hi will appear exactly 
m — 2 times in the differential equation. 

Later we will explain how to use this differential equation to obtain a numerical 
solution to Problem O 

Goluzin gave a more general result where the problem is to maximize I/'-"'' (0)| for 
any n > 1. An account of his work is in [Gol69[ Chap. 4]. By using this description 
and after considerable work, Kuzmina in |Kuz82| computes the extremal domain in 
the case of three points and in [Fed84| this is extended to four points with a certain 
symmetry (two of the points must be symmetric with respect to a line that passes 
by the other two). 

Later on, Tamrazov found an explicit solution for the problem of n points. The 
general solution to Problem [2] is, according to [Tam05j . of the form 



a<£V\{l} veW{r) I3£W^ 



l/r(v) 



where a and P are a finite number of points in T and t(v) is a positive integer. This 
result does not seem to be completely clear, because the function / corresponds to 
a Schwartz-Ghristoffel formula, thus /(T) is going to be a collection of straight 
segments (one of them going to 00) but, even when we have only three points, in 
most cases (except of very symmetric ones) the solution to the Polya-Tchebotarov 
problem are not straight lines. 
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Nevertheless the main idea of Tamrazov paper, that all solutions can be exactly 
parametrized by a finite planar graph is indeed correct. We will give a different 
proof of this fact. Our approach although it will not yield an "explicit" formula 
it will be constructive and it is possible to implement a numerical algorithm that 
produces an approximation to the solution of the Polya-Tchebotarov problem. 

On Sections [2] and [3] we prove that all solutions arc codified by a "nested par- 
tition" which arc defined there. This is more convenient for us, although it is 
completely equivalent to a parametrization by graphs. To each set of points the 
extremal continuum is in correspondence with a unique "nested partition" of T and 
conversely, each "nested partition" provides a solution. Thus all the combinatorial 
data of the solution is codified in these partitions. 

Once we have a parametrization of all possible solutions we introduce in |4] a 
numerical algorithm to compute numerically the solutions (i.e. to determine the 
parameters) to the Polya-Tchebotarov problem. Wc illustrate the method making 
it explicit in the case of 3 points and 6 points (with a certain symmetry) . This last 
case is particularly interesting because it will be of use for the applications that we 
had in mind which are developed in Section [H 

The Polya-Tchcbotarov problem is rather basic, thus it is not surprising that it 
arises in connection with many other problems. The most evident case is in the 
estimates of the univalent Bloch-Landau constant, the precise formulation of the 
problem is in Section [5l This has been exploited in |COC08| where this constant 
was improved. This work is its natural continuation. Here we will provide more 
sophisticated examples and we will use the same type of domains to improve the 
estimates of two other extremal problems that were introduced in [BC94| : the 
expected lifetime of the Brownian motion in a domain with inner radius one and 
the estimation of the principal frequency of such domains. For the precise definitions 
and results, see again Section [5l 

There are other potential applications of the Polya-Tchebotarov problem which 
could benefit from our (numerical) solution. For instance the best estimates in the 
Smale mean value conjecture obtained in |Cra07| rely on the computation of the 
solution of the problem with three points. We have not pursued improvements on 
this problem. 

Acknowledgment: We are indebted to Alex Haro for illuminating conversa- 
tions about the numerical implementation of the algorithm. 

2. The parametrization of all solutions 

In view of Laurenticv and Goluzin results ,thc continuum that is a solution of 
Problem [T] form a finite planar tree with cndpoints in the points {zi, . . . , z„}. The 
remaining nodes are of order at least three. This motivates the following definition 
introduced in [Tam05| : 

Definition 1. A graph P is a Polya-Tchebotarov graph (PT-graph for short) if it 
is a finite planar tree with the properties: 

(1) All sides of the graph are linear segments. 

(2) There are no nodes of order 2. 

(3) The sum of the length of all sides is exactly 1/2. 

Wc say that F is normalized if we mark one of the vertex of the tree (a node 
of order one). Two PT-graphs arc equivalent if there is an isomctry from one to 
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the other that extends to an orientable homeomorphism of the plane. If they are 
normahzed we also require that the isometry sends the marked vertex from the first 
graph to the marked vertex of the second graph. Wc will talk of a PT-graph T to 
denote the whole equivalence class. 

Since Problem[l]is invariant by translation we can be normalize the data {zi , . . . , Zn} 
to assume that zi = 0. 

The main result is that there is a natural way of parametrizing the solutions to 
the Polya-Tchebotarov problem by normalized PT-graphs. That is for any graph 
there is associated a unique continuum that solves a Polya-Tchebotarov problem 
and conversely all solutions arise in this way. 

Let us describe how to associate a solution to each graph. First we need one 
further definition. 

Definition 2. A partition 11 of the unit circle in a finite number of intervals is 
called a properly nested partition if the following properties hold: 

(1) The intervals come in pairs of equal length, i.e. 

n = {/i,...,/„}U{Ji,..., J„} 

and I J,; I = |/,;| for all i = 1, . . . , n. 

(2) There are no nested pairs of intervals, i.e. a couple li, Ji never separates 
another couple Ik , Jk ■ See Figure [1] 




Figure 1. Not nested and nested pairs 



Since all pairs arc not nested wc can be sure that there exist at least one pair 
li, Ji of adjacent intervals. Two partitions are equivalent if one rotation sends one 
to the other. They are normalized if we mark one of the adjacent pairs of intervals. 

It is easy to see that any PT-graph provides a nested partition and conversely. 
We start from a vertex of the graph and travel through its edges directwise. To 
each edge of the tree we associate an interval in the circle of the same length (on 
the unit circle we consider the normalized length). Each edge of the tree is visited 
twice, once for every side. We consider the pairs of intervals /i, Ji to intervals that 
correspond to different sides of the same edge of the tree. 

Finally for any given properly nested partition there is associated an involution 
r defined on the circle (except in a finite number of points corresponding to the 
end points of the intervals). Two points x^y are related by the involution r if a; 
belongs to the interval li and y belongs to its pair J,. The definition of r in each 
of the intervals is the reflection along the diameter of the disk that passes halfway 
in between the pair of intervals as in Figure [2] 

We are going to prove a "welding" type theorem: 
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Figure 2 . The involution r restricted to a pair of intervals 

Theorem 3. For any given properly nested partition II of T and its associated 
involution r there is a conformal map fjj : B> C\T where T is a finite union 
of analytic arcs T = U"^]^7i such that the image of any pair 1^, d Yl is one arc 
jk ifid fnix) = fu(T{x)) for all x in the intervals. Moreover /n is unique up to 
postcomposition with automorphisms of C 

Remark 2. If we compose /n with a translation we obtain a function that satisfies 
the Laurentiev conditions of Theorem [TJ Thus for any 11 we get a solution to 
Problem [21 The converse is even more clear. Laurentiev theorem shows that the 
boundary of the extremal domain dVl is a tree, that is homeomorphic to a rectilinear 
PT-graph F. The length of each edge of F is one half of the harmonic measure of 
the corresponding edge of dfl. The conformal map g : Q, ^ D with g{fS) =0 gives 
a partition 11 of the unit circle. By property ([3]) of Laurentiev theorem, we get 
g-^ = fn- 

3. Proof of the welding Theorem 

Given the partition wc will proceed to construct the mapping in a finite number 
of steps. In each step the following lemma is the key 

Lemma 1. Given two adjacent intervals I, J d T in the circle, such that T 7^ 
{/ U J} and a quasisymmetric homeomorphism t : I J that fixes the common 
point p, there is a simple Jordan arc 7 C B with one endpoint at 1 and a conformal 
mapping / : B ^ B \ 7 such that /(O) = and f{x) = f{T{x)) G 7 for all x G I. 
The mapping f and the curve 7 are unique, they depend on I, J and t. 

Proof Let /q {e**; t £ [0, 7r/2]} and Jo = {e**; t G [-7r/2, 0]}, and let tq : /o Jo 
be defined as tq{z) — z. If these were the data of the problem, it will be readily 
solved by the mapping fo{z) — ... that maps D to B\ [r, 1). In the general situation, 
there exists an homeomorphism of the circle T : T ^ T such that T{I) = Jq, 
T(J) = Jo and tq{T{x)) = T{t{x)) for all x G /. In /, the map T is defined 
linearly. On J we use as definition T{y) = To{T{T~^{y))) and outside / and J we 
define it linearly. Since r is asymmetric then T is quasisymmetric. In general T 
has the same regularity as r. By the Beurling-Ahlfors extension theorem |BA56| 
there is a quasiconformal homeomorphism of the disk that extends T. We will still 
denote it by T. 

Let g = fo°T- The map g is mapping / and J to the arc [r, 1] in such a way that 
g{x) = gijix)) for all x £ I. It is not a conformal, but it is a quasiconformal map 
because T is quasiconformal and /o is conformal. This can be corrected by solving 
a Laplacc-Beltrami equation. We want an homeomorphism \I> : B ^ D such that 
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fc = '5 o g is conformal. Thus fcj = 'I'jgl + '^zQz = 0. There is always a solution 
to this Laplacc-Beltrami equation that is an homcomorphism from D to B) since 
\gz/gz\ < k < 1 almost everywhere by the measurable Riemann mapping theorem. 
Finally we compose k with an automorphism of the disk k and the desired function 
is f ^ K o k. The automorphism k is chosen to make sure that /(O) = and the 
endpoint of /(/) = 5" ([r, 1]) is at 1. The regularity of / at the boundary is as good 
as that of 'i> and that itself is determined by the regularity of r. 

If there were two curves 70 and 71 and two maps /o and /2 with the same 
property then g = fio f^^ is a conformal mapping from E* \ 70 to D \ 71 such that 
5(0) = 0. Moreover g extends continuously to 70 because for any point in 70 the 
preimage by /q are two points {x, y\ in T that are related by r, thus /i(x) — fi{y), 
therefore g extends continuously to 70 and therefore it extends analytically, thus 
g{z) = e^^z. Since both 70 and 71 start have an endpoint in 1, then g{z) = z and 
fo - /i. □ 



With the same proof, mutatis mutandi, we deal with the case T = / U J and we 
obtain 

Lemma 2. Given two adjacent intervals I, J d T in the circle, such that T = 
{I U J} and a quasisymmetric homcomorphism t : I —t J that fixes the common 
points, there is a simple Jordan arc 7 C C with one endpoint at 1 and the other 
at 00 and a conformal mapping / : D ^ C \ 7 such that /(O) = 0, /'(O) > and 
f{x) = f{T{x)) G 7 for all x G /. The mapping f and the curve 7 are unique, they 
depend on I, J and r. 

Proof of the theorem. We take any pair of adjacent intervals {Ik, Jk) in the partition 
n corresponding by the involution r. There are always adjacent pairs because the 
partition is properly nested (they correspond to edges with and endpoint in a vertex 
of the graph). Applying the Lemma we find a conformal mapping /i that welds 
together the pair of intervals in a curve 71 C D. The mapping /i induces a new 
partion Hi of T, = fi{Ij) and Jj = fi{Jj) for all j except for the pair (Ik, Jk) 
that was welded together. This new partition is again correctly nested because 
the order is preserved except for a pair of adjacent intervals that "collapses" . The 
number of pairs of intervals is one less than in Hi. The intervals in each pair are 
no longer of the same size but nevertheless the map /i induces a new involution on 
them, Ti = fiOTo f[~^\ 

Now we repeat the procedure. We take any other new pair of adjacent of intervals 
of the new partition IIi corresponding by ri and we glue them together. This 
can be done by a mapping /2 applying the Lemma because ri is quasisymmetric 
(it is in fact piecewise real analytic). In this way we get again a new involution 
^2 = /2 o Ti o /2 and a new nested partition 112. 

In this way we keep gluing pairs of intervals until we are left only with two 
intervals and an involution r„ that relates them. In this last step we use Lemma [2] 
to get /„. The final conformal mapping is /n = /ji o • • • o /i- 

There is basically only one such map /n (except for composition with maps 
of the form az + b). The proof is as in Lemma [l] suppose there is another such 
map gu- Let F = /n(T), then h = gu ° fj^ is a one to one analytic mapping 
h : C\T ^ C Moreover since the preimage of any regular point in F are two 
points in the circle that are related by the involution, and g maps both points to 
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the same point, then h can be extended continuously to a conformal map from C 
to C, thus h{z) = az + b. □ 

4. Numerical Algorithm to find solutions 

As we mentioned above, in order to implement numerically an algorithm to find 
the solution to Problem [21 we used an important property of the extremal domain 
and the differential equation obtained by Laurenticv. Denote by ri„ the desired 
extremal domain for the Problem [2] in case of ri + 1 points {ai, . . . , a„, oo}. Let 
/ : D ^ 0„ be the conformal map such that /(O) = 0. We know that / satisfies 
the following differential equation 



/(-) J uUim-b 



where the parameters bj are unknown and C = j-j,,-^ {-ak) ' Using the solutions of 
this differential equation and the last property of Theorem [T] we have implemented 
the resolution of Problem [2] for some cases of n. The system becomes more deli- 
cate as n increases, the combinatorics and the dimensions of the systems to solve 
become bigger. We will show in detail the solution in the case n = 3 to illustrate 
the method and n = 6 with some extra symmetry, because this will be enough 
for the applications that we have in mind. The code where this algorithm is im- 
plemented (for four points and 6 points with symmetry) can be downloaded from 
_http://www. maia.ub.es/cag/code/tchebotarev/ 




Figure 3. Sketch of the extremal compact for three points 

4.1. Case of 3 points. Let's start with 3 points. Assume that we have 01,02,03 
three points such that 03 = 00 and 01,02 7^ 0. Without loss of generahty we will 
always assume that /(I) = 00. In the case of three points the extremal domain is 
very clear (see Figure [3]) . We only have one unknown parameter denote it by b in 
the differential equation ^ that reduces to: 

^2) j'(^)2 ^ (jili^) - ai)(/(2) - 02) f{zf 



f{z)-b z2 

Recall that to any regular arc afi of 9f2 there corresponds two arcs with equal 
lengths on the unit circle. So this means that we have the configuration on the unit 
circle shown in Figure [H where /(O) = 0, the arcs (l/3i), (/Jsl) are mapped to the 
arc 00b, Piai and ai/32 are mapped to the arc 601 and the arcs /32Q!2, a2/33 into 602. 
Note that f ie'"') = Oi for i = 1,2 and /(e**) = 6 for i = 1,2,3. 
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Figure 4. Configuration for n = 3 



The solution of the problem can be viewed as the solution of a system of non- 
linear equations. If we know the value of /'(O) and b we can compute the coefficients 
of the mapping / using the differential equation 0. We computed also the values 
of Ui.Pi. Note that as the arcs l/3i,/33l must have same length, we get /Js = 27r — /3i 
(we will always take the angles in the range [0, 27r)). So, we have 6 real unknown 
parameters in our problem: Re(/'(0)), Im(/'(0)), Re(5), Im(6), /?2 and using the 
last property of Theorem [1] we can impose the following three complex equations 

j^^i(a2+P2)/2-^ ^ y(gi(a2+/33)/2)^ 

We used a the hybrid method to find an approximation of the roots of the system 
(see jPowTOj for more details of the method) . 

To apply the root-finding method we need to evaluate f{e^'^) for any 7 S [0, 2tt) \ 
{ai,a2, P2, 03}- For that, denote z{t) = f{te^''). We know that z(0) = and 
z'(0) = /'(0)e*". Note that z(l) = /(e^"). We can get the difTerential equation 
satisfied by z{t) and solve it to obtain the value in time t = 1. We get 

{z{t) - ai)iz{t) - a2) z{t)^ 



(3) 



z'(t)2 = 



{z{t)-b) 

Note that this equation only defines z' [t) up to a sign, we will deal with this problem 
by analytic continuation. Once we fix the derivative at the origin there is a single 
analytic branch that solves the equation. To solve it we used the Taylor integration 
method which allows us to integrate the singularity in t = Q. As / is conformal, 
we know that z{t) = zit + Z2i^ + . . ., where zi = /'(O)e*'''. Now if we do the 
calculations in the equation ([3]) we get a recurrence for the coefficients till the order 
we want. Now we can estimate the radius of convergence of the obtained series. 
And therefore proceed using Taylor method to integrate the differential equation 
until t = 1. Hence we will be able to impose the equations to solve our problem. 

4.2. Case of 6 points with symmetry. Now consider the case of 6 points 
ai,a2,...,a6 such that ag = 00, 03 e M and 05 = 01,04 = d2. The extremal 
compact in this case may be of two types (see Figures [H] and [7]) . So, we have two 
type of configurations in this special case 

(1) 13\ ai Pi a2 /3| as /3| ^4 /3| Pi as Pi 2tt 

(2) Pi Pi ai Pi a2 Pi Pi as Pi Pf Pi as Pi Pi 2^, 

where /(e*'^^) = bj for j = 1,2 and k = 1,2,3,4, /(e'"-) = o^ for i = 1,2,3,4,5 
and /(I) =00. Using the symmetry we can do some reductions to get a system of 
equation with less dimension. For example, we can always assume that the point 
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Figure 5. Structure of the extremal domains for n = 6 with symmetry 

a3 G M, so that 61,62 E M for the first configuration and 61 e M for the second 
one. Note that in this last configuration, by the symmetry of the problem, 63 = 62- 
Moreover as this is a symmetric case, /'(O) must be real and ~ tt. 

First configuration. In this case we have 7 real unknown parameters: Re(/'(0)), 
Re(6i), Re(62), Pi, [3\, f3f, Pf and we can impose the following equations 

'lm(/(e*^'iV2.0))^0 

< Im(/(e*('32+'3?)/2)) = 

^Im(/(e*("3+'32)/2)) ^ 

Second configuration. We have 8 real unknown parameters: Re(/'(0)), Re(6i), 
Re(62), 1111(62), Pf, Pi^Pi impose the following equations 

'lm(/(e*'5j/2.O))^0 

J(gi(ai+/3?)/2-) ^ y^gt(Qi+/3|)/2-) 

< J(gi(Q2+/32)/2) = y(gt(Q2+/9^)/2-) 

Im(/(e*("^+'33)/2)) = 

In Figures [6] and [7] we show an extremal domain for some 6 points for each con- 
figuration. As in the last case, these figures represent a conformal map g from the 
complement of B onto such that g{oo) = 00. 

Remark 3. Note that the solution of the problem depend continuously on the pa- 
rameters ai, so if we have one solution for some given set points we can do continu- 
ation to reach to any other set of points (with the same topological configuration). 
This has been used and we did the classic continuation (i.e. for the new set of 
points we take as a initial condition the solution of the last set of points). 

Remark 4. In the implementation of the method, we found a problem when the 
distance between the arcs on the unit circle is very small, we can't integrate properly 
the differential equation because we are near the poles bi. However this can be 
overcome by a change of variables. This has been implemented for the special 
case of the first configuration mentioned in the case of 6 points. In fact, for the 
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Figure 6. Extremal domain for n = 6 with symmetry (configu- 
ration 1) 
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Figure 7. Extremal domain for 71 = 6 with symmetry (configu- 
ration 2) 



application of the Polya-Chcbotarev we only needed the values of /'(O), bi,b2- This 
data is enough to obtain the expansion in series of the mapping /. This can be done 
in the following way: given the points a^, we have an initial guess for the unknown 
values. So, we compute the critical orbit starting at point ai till the point with 
imaginary part equal to Im(ai/2.0) and same for the critical orbit starting at bi. 
The real part of the two points obtained should match if 61,62 are the desired 
solution. So this is one real equation. Same can be done for the couple of points 02 
and 62- So we have two real equations. The last equation can be Im(/(e*°'^)) = 
(this is valid if the points ai, 02 are a bit far from the point og). 
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5. Applications of the Polya-Tchebotarov problem 

The fundamental frequency of a domain. In 1965, Endrc Makai (see [Mak65| ) 
proved the following theorem solving a problem in the study of vibrating mem- 
branes raised by Polya and Szego in their book |PS51| . In 1978, Hayman ( |Hay78| ) 
unaware of it, reproved the same result. 

Theorem 4. Let D be a simply connected domain in the complex plane. Let Ru 
be the inradius of D, that is, the radius of the largest disc contained in D and let 
Xd be the first Dirichlet eigenvalue for the Laplacian in D. There is a universal 
constant a such that 

(4) A,, > 

There have been many efforts to find the best constant a and to identify the 
extremal domain for a. Makai's proof also shows that the best a satisfies 1/4 < 
a < 77^/4. 

The following lemma is useful for giving upper bounds for this constant (see 
jBC94( Lemma 1.2] for a proof): 

Lemma 3. Let Jo{r) be the first Bessel function and jo the smallest positive zero 
of Jq. Assume that D is a simply connected domain. Then 

\d < jo inf ' 

wh 



!',4{3ory--'dr 



/o J§{jor)rdr 

and the infimum is taken over all conformal mappings F{z) = Yl^=o o-n^" from the 
unit disc onto D. 

In |BC94| Baiiuelos and Carroll proved that 0.619 < a < 2.13 and provided 
examples of domains which are close to the extremal domain. They did this relating 
this problem to two other extremal problems: 

The expected lifetime of a Brownian motion. Let Bt be the Brownian motion in D. 
Let TD ~ inf {t > Q : Bt ^ D} be the first exit time of Bt from D. Let us denote by 
Ez {jD ) the expectation of td under the measure of the Brownian starting at the 
point z in D. It is known that there is a universal constant b such that, whenever 
_D is a planar simply connected domain, 

(5) sn^i Ez{td) <bRl. 

z&D 

As before, we want to know the best value of h and the extremal domain for this 
last inequality. It is a fact that if D denotes the unit disc of radius Rd then 
sup^g£,£'^(r£,) = i?|,/2. It is known that 1.584 <b< 3.228 (see jBC94| ). In order 
to give an improved lower bound for b it is useful to know the following result (see 
[BC94i Lemma 1.1.] for further details): 

Lemma 4. Suppose that F{z) = X)^o ^nz" is a conformal mapping from the unit 
disc onto a simply connected domain D with F(0) ~ zq. Then 



oo 



71=1 
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The univalent Bloch-Landau constant. If / is an analytic and one to one mapping 
from the unit disc, then there exists a universal constant U such that 

(6) %(D) >Z^|/'(0)|. 

This means that the image of the unit disc under any conformal map / contains 
discs of radius less that U\f'{Q)\. Note that from Koebe's 1/4-theorem, we know 
that U > 1/4. The best value of U is known as the univalent or schlicht Bloch- 
Landau constant. We can reformulate this problem in terms of the density of the 
hyperbolic metric. If / is a conformal mapping from the unit disc such that /(O) = z 
then the density of the hyperbolic metric is cr{z; D) = 1/|/'(0)|. So we have the 
following inequality 

(7) ao:^ mla{z;D)>^. 

where c := U. From many years there have been efforts to find bounds for U ~ c. 
This constant was introduced in 1929 by Landau [Lan29| . who proved that U > 
0.566. Reich improved this bound in [Rei56| {U > 0.569) and Jenkins in [Jen61| 
gave U > 0.57088. Many other gave some improved bounds. There are many 
domains proposed as the candidate for the extremal domain in order to obtain 
upper bounds for the Bloch-Landau constant. For example, Robinson in |Rob35| 
proved that U < 0.658, Goodman in |Goo45j that U < 0.65647 and in |BH85j Beller 
and Hummel proved that U < 0.6564155. Finally in [COC08| this bound has been 
improved to U < 0.6563937. In this last result the resulting domain had all the 
inner boundary harmonic symmetric with respect to the origin, see Figure [51 We 




Figure 8. Structure of the domain in [COC08| 



will use a slight modification of this last domain, which is still harmonic symmetric, 
to give an improved bound of the Bloch-Landau constant. 

As a normalization, we will take domains D with inradius 1. We will give 
improved bounds for the constants appearing in the three problems explained above. 

Baiiuelos and Carroll in |BC94| conjectured that the extremal domain is the 
same for all the three problems. When we restrict these problems to the class of 
convex domains this is true. In our work, we will see that a similar domain improves 
the bounds for all the three problems. 
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5.1. Bloch-Landau constant. If /i is a mapping of the complement of a compact 
set E onto the complement of the closed unit disc, it can be expanded (up to a 
rotation) as 

M.) = ^ + 0(l),z^oo. 

So we can relate the problem of the extremal domain for the Bloch-Landau 
constant with Problem [2] in the case of 6 points because minimizing the capacity is 
equivalent to increasing the derivative at the origin. It is known (see jCarOSj ) that 
the arcs making up the extremal configuration must be harmonically symmetric 
at infinity. We will work with domains il = ^zi,z2,R ^ in Figure [H] where R is 
bigger than 4 and we chose the arcs 71 , 72 , 7i , 72 so that this domain is harmonically 
symmetric with respect to 0. If 5 is a conformal map of B onto ^zi,z2.R with g(0) = 




Figure 9. Domain ilzuZ2,R 

then f{z) = Zy/g(z^)/z^ is a conformal map of B onto the domain Dw-^^^^n shown 
in Figure 1101 The arcs in this last domain are harmonically symmetric. We can 
compute the derivative of / at 0: |/'(0)| = Aj/|g'(0)|. We need this domain to have 
inradius 1. Later on we will explain the construction of the domain and the way 
to get inradius 1. Let k(z) = z/{\ — z)^ be the Kocbe mapping from the unit disc 
onto the complex plane slit along the negative real axis from minus infinity to -1/4. 

Proposition 1. Let f be the conformal map o/B onto -Dtui,iu2,i? such that /(O) = 0. 
Then taking z\ = w\ and Z2 = , 

^ = ^m-8)-V'(l)|cap(i?) 

where ipi^) = —l/k{z/R'^) and E is the continuum with minimal capacity contain- 
ing 6 given points with symmetry (see FiQure \TJ\ ) . 

Proof. Using the above notations, consider the map 

_j^{z)-jm_ 

- |^(-8)-^(l)| 
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ai 



Figure 1 1 . The continuum E 



which maps ^lzi,z2,R onto the complement of the continuum E with ai — (f>{zi) and 
02 = (f>{z2) (</>(!) = 0, (?!)(— 8) — 1,0(0) = oo). Note that the harmonic symmetry 
of the arcs are preserved since each mapping can be extended continuously to all 
internal boundary arcs and the domains involved are symmetric with respect to the 
real axis. Let h be the mapping of the complement of E onto the complement of 
the unit disc. Then we can define the map G of Qzi,z2,R onto D as 

G{z) ^ 



We can calculate the derivative G'(0) by computing the power series of G: 

R^h'ioo) " • 

The capacity of E can be computed numerically as we explained in 14.21 Now 
F /^^ = Z'^G{z^)/ and the proposition is proved. □ 
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5.2. Construction of the domain Dwi,w2,R- Now we will explain how to obtain 
the desired domain in order to have inradius 1 and tlie derivative |./'(0)| as big as 
possible. In what follows, Dp denotes a disc of radius 1 centered at p and Cp — dDp. 

The domain is constructed in stages. Let -D(0, R) be a disc centered at the origin 
with radius R> 1. First we remove from this disc three radial slits that start from 
the cube roots of the unity. Then, we remove three further slits starting at two 
times the cube roots of -1 (these are the first two stages of Goodman's domain). 
Now, as we need inradius 1 we need to put some point in order to not to have 
discs of radius bigger than 1. Let Pi ~ (1 + V 2-^3 — 3, 1) and denote by Ci the 
circle centered at Pi with radius 1. We have to put some point in Ci so that this 
circle can't increase. Let wi be a point in this circle. Denote by C2 the circle of 
radius 1 tangent to the halfline of argument tt/B containing the point wi. Let C3 
be the circle of radius 1, tangent to |z| = i? and to the halfline of argument 7r/3 and 
denote by the intersection point of Ci and C2, P2 and P3 the centers of C2, C3, 
respectively (see Figure [T2|) . Now let Fi and F2 be the curves at distance one of 




Figure 12. Election of wi and UI2 

71 and 72 (i.e. given a point ^i{t) £ 7^, let Vn be the normalized orthogonal vector 
to 7i(t), then the corresponding point in the curve Fi is Ti{t) = ji{t) ± Vn- Let q 
denote the intersection of Fi and F2. One sufRcient condition to have inradius one 
is \q\ > R — 1. The idea to prove this is to cover all the points {x,y) so that they 
can't be centers of circles (contained in Dw^^w2,r) with radius bigger than one. Let's 
show it when our points are located in the sector between the segment and 
the halfline of argument 7r/3. If z = (x, y) is a point such that y < l,\z\ > R — 1 or 
dist(2:, T) < 1 (where T is the halfline of argument 7r/3) then obviously we cant have 
a circle centered at such point with radius bigger than one. Let and Dy^^ be the 
discs of radius one centered at the points wi and 'W2, respectively. If z G D^i UD^^^ 
then Dz contains one of the points wi or W2- Therefore, the region D^,^ U D^^^ is 
covered. The only risky region is the one between dD^i, dDu,2, l^j = i? — 1 and 
{z = (x,y)\y ~ 1}. But this space is covered by the region delimited by the curves 
Fi U F2 due to the hypothesis (see Figure [13]). 
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Figure 13. Prohibited zones 

5.3. Results. We have computed the bounds for aU three problems explained be- 
fore. To construct the point wi wc move on the real axis x >= 1 + \/ 2-\/3 — 3 and 
define the point P2 and then wi. After this all gets determined. So given x, first we 
find the biggest R such that |(?| > 1 — R, (because the derivative at the origin will 
increase with the radius R) and after that we compute the bounds of the constants 
explained in the three problems. The results obtained are: 

(1) For the Bloch-Landau constant, the best upper bound has been found for 
X = 2.1383799965243 and R = 5.1195152501 and the improved bound is 

U < 0.656319277272. 

The domain obtained is shown in Figure 1141 

(2) Computing the coefficients of the conformal mapping obtained for the do- 
mains Dwi.w2.R, tlic improved upper bound for the fundamental frequency 
has been found for x = 2.1282995811037759 and R = 5.10223601895443 
and it is 

a < 2.0907934752309. 

(3) The improved lower bound for the expected life time of a Brownian motion 
has been found for x = 2.174447128952 and R =■ 5.1836816989 and it is 

h > 1.670724582110. 

In all the computations the error estimates that we required are of the order of 
10~^^, and we are pretty confident on the correctness of the 10 first digits on the 
bound, but we have not done a rigorous error analysis. 
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